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The  advent  of  supermaneuverable  flight  trajectories  for  high  performance  aircraft 
opens  up  new  research  areas  involving  the  human  factors  effects  on  pilots  when 
they  fly  these  unusual  flight  scenarios.  There  exists  a  need  to  conduct  studies  in 
motion  simulators  such  as  centrifuges  to  investigate  how  these  unusual  motion  fields 
effect  humans  both  physiologically  and  how  they  impact  the  ability  of  a  pilot  to 
perform  a  mission.  The  dynamic  control  of  a  centrifuge  simulator  to  emulate  super¬ 
maneuverable  flight  trajectories  is  studied  in  this  paper  by  considering  the  centri¬ 
fuge  simulator  as  a  3  revolute  axes  robot  manipulator. 
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A  Study  of  Supermaneuverable 
Flight  Trajectories  Through  Motion 
Field  Simulation  of  a  Centrifuge 
Simulator 

The  advent  of  supermaneuverable  flight  trajectories  for  high  performance  aircraft 
opens  up  new  research  areas  involving  the  human  factors  effects  on  pilots  when 
they  fly  these  unusual  flight  scenarios.  There  exists  a  need  to  conduct  studies  in 
motion  simulators  such  as  centrifuges  to  investigate  how  these  unusual  motion  fields 
effect  humans  both  physiologically  and  how  they  impact  the  ability  of  a  pilot  to 
perform  a  mission.  The  dynamic  control  of  a  centrifuge  simulator  to  emulate  su¬ 
permaneuverable  flight  trajectories  is  studied  in  this  paper  by  considering  the  cen¬ 
trifuge  simulator  as  a  3  revolute  axes  robot  manipulator. 


Introduction 

Supernianeuverability  (Clark,  1987;  Cook,  1989;  Cord,  1990; 
Hague,  1981;  Schefter,  1989;  Sweetman,  1990)  refers  to  the 
unusual  flight  trajectories  presently  being  investigated  with 
high  performance  fighter  aircraft.  Present  fighters  can  be  made 
to  fly  these  unusual  flight  maneuvers  by  several  known  meth¬ 
ods;  using  a  high  angle  of  attack  (a  >  SO  degrees),  thrust 
vectoring,  vortex  manipulation  for  control,  or  possibly  using 
differential  deflection  of  stabilizers  or  canards.  Most  super- 
maneuverable  flight  trajectories  have  not  actually  been  flown 
in  aircraft  but  studies  with  1/7  th.  size  prototypes  are  presently 
being  conducted  to  investigate  structural  problems  and  the 
measurement  of  the  linear  and  rotational  accelerations  and 
forces  that  would  act  on  a  pilot  (Gal-Or,  1990)  during  a  spec¬ 
ified  maneuver.  In  his  work  a  set  of  all  pitch  and  roll  angles 
of  a  specified  aircraft  (F-15)  were  determined  where  both  con¬ 
trollable  and  stable  flight  will  occur.  This  region  of  aircraft 
attitude  is  clearly  a  function  of  both  the  altitude  and  speed 
due  to  the  structural  problems  and  controllability  issues  as¬ 
sociated  with  flying  an  aircraft  in  this  unusual  regime. 

One  particularly  well-known  supermaneuver  is  the  “Cobra 
Maneuver”  which  has  been  recently  flosvn  at  airshows  (Avia¬ 
tion  Week,  1990)  by  both  Soviet  and  U.S.  high  performance 
aircraft.  Recent  work  (Chiang  et  al.,  1990)  has  addressed  the 
problem  of  controller  design  in  these  unusual  flight  scenarios 
in  which  there  exists  a  compromise  of  high  angle  of  attack  to 
controllability  and  stability.  In  their  study  an  FA/ 18  was  uti¬ 
lized  because  of  its  use  of  thrust  vectoring  as  a  procedure  to 
produce  the  supermaneuver.  Their  motivation  was  based  on  a 
previous  study  (Kalviste,  1986)  in  which  it  was  shown  that  a 
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pilot  must  control  both  the  standard  P  and  Q  angles  simul¬ 
taneously  to  change  the  aircraft’s  body  axis  roll  angle.  This  is 
because  neither  the  body  Af-axis  nor  the  velocity  vector  (wing 
-T-axis)  can  be  used  as  the  roll-axis  reference.  Both  these  vari¬ 
ables  change  direaions,  themselves,  during  the  supermaneu¬ 
ver. 

One  main  issue  which  has  not  been  addressed  is  how  such 
unusual  flight  scenarios  will  affect  pilots  in  their  ability  to 
continue  to  perform  tasks  during  the  mission.  To  understand 
better  these  human  factors  issues,  it  is  necessary  to  conduct 
research  with  pilots  on  a  ground  based  motion  simulator  that 
will  reasonably  replicate  the  motion  field  experienced  by  pilots 
flying  these  unusual  maneuvers.  The  term  “reasonably  repli¬ 
cate”  is  mathematically  formulated  in  this  paper  using  a 
Damped  Least  Squares  numerical  procedure.  A  motion  sim¬ 
ulator  (large  centrifuge)  exists  at  Wright  Patterson  Air  Force 
Base,  Dayton,  Ohio  in  the  form  of  a  large  three-degree-of- 
freedom  device  which  is  used  to  test  pilots  and  equipment.  The 
analysis  of  its  control  to  produce  trajectories  to  emulate  su- 
pemaneuverablc  aircraft  will  be  conducted  in  this  study  by 
analyzing  the  centrifuge  simulator  within  the  context  of  a  large 
robot. 

A  centrifuge  simulator  differs  from  traditional  robots  be¬ 
cause  of  three  fundamental  reasons.  First,  its  size  is  prohibi¬ 
tively  large  (190  tons  in  weight  and  25  feet  in  height)  but  it 
can  be  considered  a  robot  with  3  revolute  joints  (all  axes  have 
360  degrees  of  rotation).  Second,  the  centrifuge  has  a  simpli¬ 
fication  (due  to  balancing  conditions)  of  having  the  same  dy¬ 
namic  payload  at  all  points  during  its  dynamic  operation  (the 
inertia-mass  matrix  is  constant  for  all  time)  This  is  a  great  aid 
in  simplifying  the  computations.  The  third  way  the  centrifuge 
simulator  differs  from  traditional  robots  is  due  to  the  fact  that 
its  goal  is  to  create  an  unusually  different  motion  field  at  a 
panicular  point  on  the  end  effeaor  (at  a  point  on  a  seat  within 
the  centrifuge).  This  motion  field  will  consist  of  a  specified 
linear  acceleration  (a  three-dimensional  vector)  which  changes 
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as  a  function  of  time  and  also  a  particular  rotational  motion 
field  of  a  fixed  coordinate  system  at  the  rigid  body  defined  at 
the  endpoint  on  the  seat  of  the  centrifuge  where  the  pilot  is 
located.  This  rotational  motion  field  may  also  vary  as  a  func¬ 
tion  of  time.  The  study  of  this  type  of  problem  differs  from 
the  traditional  control  theory  approaches  in  robotics  where  it 
is  desired  to  control  the  position  of  the  end  effector  or  the 
force  at  the  end  effector  in  contact  with  the  environment.  First 
it  is  necessary  to  formulate  the  problem  of  interest  as  it  applies 
to  a  centrifuge  simulator. 

Formulation  of  the  Problem  of  Interest 

Let  A”  be  a  3  X  1  position  vector  in  Cartesian  space  and  d  a 
3  X  1  joint  space  vector.  The  forward  kinematics  relationship 
between  X  and  6  can  be  specified  via: 

X=-f{e)  (1) 

where  the  function  /(.)  (cf.  Repperger,  1990)  for  a  centrifuge 
motion  simulator  is  an  extremely  complicated  function  of  the 
3  joint  variables  (/  =  i,  3).  If  A"  is  a  linear  position  vector 
and  denoting  its  3x  1  dimensional  velocity  vector  as  F,  then; 

V=^X  =  Jld  (2) 

as  the  linear  velocity  relationship  and  the  dot  indicates  time 
differentiation.  Jl  is  the  3  x  3  Jacobian  matrix  associated  with 
linear  velocities. 

If  the  vector  fi  represents  the  angular  velocity  of  a  coordinate 
frame  located  at  the  end  effector  within  the  simulator,  then 
the  representation; 

Q  =  J2d  (3) 


is  derived  in  detail  in  the  sequel.  Here  J2  is  another  3x3 
Jacobian  matrix  associated  with  rotational  velocities.  Com¬ 
bining  equations  (2,3),  it  is  seen  that; 

col  [F,n]  =  coi(ji,  y2i@  =  y«  (4) 

where  y  is  defined  as  the  overall  6x3  Jacobian  for  this  system, 
yean  be  calculated  in  a  variety  of  ways  (see  Orin  and  Schrader, 
1984)  and  in  the  sequel  one  method  to  calculate  this  matrix  is 
presented  for  the  centrifuge  from  the  forward  kinematic  equa¬ 
tions. 

What  is  of  interest  in  .4ir  Force  applications,  however,  is 
the  linear  accelerations  acting  on  a  pilot  located  at  the  endpoint 
of  the  centrifuge.  This  motion  field  produces  forces  acting  on 
a  pilot  which  physiologically  stress  him.  In  addition,  it  is  nec¬ 
essary  to  separate  out  the  Coriolis  accelerations  from  the  cen 
irifugal  and  gravity  terms.  The  Coriolis  accelerations  are 
extremely  undesirable  in  both  linear  acceleration  stress  studies 
and  in  rotational  studies  involving  spatial  disorientation.  How¬ 
ever,  it  is  necessary  to  keep  careful  account  of  these  Coriolis 
terms  due  to  efforts  to  dynamically  control  out  these  particular 
terms.  To  develop  the  linear  acceleration  terms,  the  chain  rule 
is  used  in  differentiating  Eq.  (2)  to  yield: 

a=l>=yi6+yie  (5) 

where  o  is  the  3x  1  acceleration  vector  and  / 1  would  be  a 
third  order  tensor  in  this  notation.  The  term  6  is  used  in  this 
context  to  denote  a  vector  of  joint  accelerations.  Thus  to  de¬ 
velop  the  motion  field  desired  within  a  centrifuge  simulator 
for  studies  involving  man-machine  systems,  one  is  particularly 
interest  in  the  relationships  (3)  and  (5). 

What  complicates  the  problem  is  that  the  3  x  1  vector  a  on 


Nomenclature 


Of  = 

X  = 
«,(/=!  =3)  = 


e,  = 

At  = 
V  = 

Ji(,=  l  =  ,2)  = 

n  = 

6,  =  u),  = 

a  = 
C  = 


L  = 


Ri  (/=1,2)  = 
X  = 


(nr  ((=1.2)  = 


P,  (f  =  0,5)  = 


the  angle  of  attack  of  an  aircraft 
(pitch  angle) 

the  cartesian  position  vector 
the  joint  angles  of  the  centrifuge 
or  robot 

stepwise  time  integration  of  the 
joint  vector  9 
time  step 

linear  velocity  vector 
Jacobians  of  the  robot  or  centri¬ 
fuge 

angular  rotation  vector 
specific  joint  rotation  velocities 
for  joint  /  (f=  1,3) 
linear  acceleration  vector 
a  scalar  constraint  on  the  joint 
velocities  as  a  result  of  a  Euclid¬ 
ean  norm  operation 
cost  functional  to  weigh  both  ro¬ 
tational  motion  artifact  and  limits 
on  the  finite  joint  velocities  al¬ 
lowed 

the  weighting  matrices  in  the  L 
cost  functional 

adjoint  variable  of  the  cost  func¬ 
tional  L  used  to  weigh  joint  ve¬ 
locities 

residual  terms  of  the  cost  func¬ 
tional.  rl  weights  motion  artifact 
in  rotational  space;  r2  weights 
limits  on  joint  velocities 
position  vectors  describing  the  se¬ 
lection  of  coordinate  frames  using 
the  Denavit-Hartenberg  formula¬ 
tion 


LI,  L2,  Lx,  Ly, 

Lz,  Lt,  e„  and  8,  = 


[X.  Y„  Z,)(/  =  0,5)  = 


T1  = 


R,  Px\,Py\,Pz\  = 


(mi  jm  (m  "0i5)  — 

Aj(t),  Bj(t),  and 

Cj(t)  a  =  0,3)  = 


9x(t),9y(t),9z(t)  = 


^j{t),(j-x,y,z)  = 


fixed  quantities  of  the  centrifuge 
simulator  including  arm  and  seat 
lengths  and  seat  rotations  to  de¬ 
scribe  the  Denavit-Hartenberg  pa¬ 
rameters  in  the  transformation 
matrix 

describe  the  coordinate  frames  se¬ 
lected  in  the  Denavit-Hartenberg 
transformation 

the  transformation  matrix  from 
frame  i  to  j  using  the  Denavit- 
Hartenberg  formulation 
notation  in  Denavit-Hartenberg 
formulation  where  R  represents 
the  rotation  matrix,  and  the  Pj 
(/'  =  x'l,  >’1,  zl)  are  position  vec¬ 
tors 

unit  vectors  in  the  wth  selected 
coordinate  frame 

coefficient  matrices  of  the  unit 
vectors  in  the  various  coordinate 
frames  of  interest  to  describe  ro¬ 
tational  motion 
actual  (body  axes)  rotations  a 
pilot  sees  when  he  sits  in  the  cen¬ 
trifuge  simulator  or  in  an  actual 
aircraft  during  the  supermaneuver 
actual  (body  axes)  rotational  ve¬ 
locities  a  pilot  sees  when  he  sits  in 
the  centrifuge  simulator  or  in  an 
actual  aircraft 


Journal  of  Dynamic  Systems,  Measurement,  and  Control 


JUNE  1992,  Vol.  114/271 


the  left-hand  side  of  (5)  is  specified  as  a  function  of  time. 
Simultaneously  the  3  x  1  vector  0(0  in  (3)  is  now  determined 
if  the  0,  variables  in  (5)  are  initially  specified  by  the  require¬ 
ments  for  a(/).  Thus  for  the  simulation  of  certain  motion  fields, 
there  exists  situations  where  the  problem  may  be  over  or  under 
determined  and  one  must  be  careful  to  determine  if  the  linear 
or  rotational  motion  field  specified  provides  a  well  posed  prob¬ 
lem  for  the  given  simulator.  It  is  necessary  to  define  two  prob¬ 
lems  of  interest. 

The  first  problem  of  interest  would  be  the  seleaion  of  the 
joint  commands  d,{t)  such  that  the  linear  acceleration  vector 
a{t)  at  the  end  effectoi  would  match  a  desired  motion  tra¬ 
jectory.  The  second  problem  would  be  the  selection  of  the  joint 
commands  6,{t)  to  match  the  rotational  time  histories  of  the 
coordinate  frames  located  at  the  end  effector  to  emulate  a 
specified  motion  field.  For  simplicity  the  second  problem  will 
be  considered  in  this  paper. 

For  the  robotics  problem  of  interest,  it  is  desired  to  find 
0  to  achieve  a  particular  motion  trajectory  X(t)  which  is 
specified  apriori.  To  relate  this  procedure  to  a  practical  prob¬ 
lem  (Repperger,  1990, 1991)  involving  motion  base  simulation 
of  a  specified  aircraft  trajectory  on  a  simulator,  we  define  as 
the  problem  of  interest  (e.g.,  using  only  Eq.  (3)  when  J2  is 
square)  to  solve  for  d,  when  Q  may  be  specified,  i.e.: 

0  =  lJ2]-'a  (6) 

If,  however,  J2  is  in  a  nearly  rank  deficient  configuration,  this 
leads  to  excessively  high  joint  velocities  0 .  There  may  exist  a 
practical  constraint  of  the  form; 

(7) 

to  make  such  a  problem  realizable.  The  norm  used  in  Eq.  (7) 
refers  to  the  standard  Euclidean  norm. 

The  Damped  Least  Squares  numerical  method  provides  an 
ideal  way  to  convert  the  constraint  problem  specified  by  (6-7) 
into  an  optimization  procedure. 


Damped  Least  Squares 

This  type  of  algorithm  is  ideal  when  the  Jacobian  may  be¬ 
come  nearly  singular  (Wampler,  1986;  Maciejewski  and  Klein, 
1988).  To  understand  how  this  procedure  handles  the  singu¬ 
larities  of  the  Jacobian,  the  objective  is  to  minimize  L  specified 
by: 

=  (8) 

where  X  >  0  is  selected  to  be  constant  in  this  paper  and  to 
minimize  joint  velocities  if  the  Jacobian  becomes  singular. 
Thus  the  fidelity  of  the  motion  field  (denoted  by  lX-J0i^) 
is  compromised  by  trading  off  a  cost  associated  with  a  weight¬ 
ing  in  X^19l^.  This  is  appropriate  for  a  centrifuge  simulator 
in  which  finite  torque  on  the  axes  limits  the  0  values.  Also, 
since  the  energy  in  the  simulator  is  proportional  to  0 ,  this  is 
consistent  with  the  finite  energy  available  to  produce  motor 
torques. 

By  minimizing  6  of  Eq.  (8),  this  yields  the  augmented  equa¬ 
tion: 

col[y,  X/]fl=col(A',  0)  (9) 

which  requires  solving; 

(J^J+X^f)0=J^X  (10) 

or 

0=(fj-^\^I)-'j^X  (11) 

as  the  inverse  kinematics  solution.  Thus  if  X  =  0  in  Eq.  (11), 
this  reduces  to  the  results  specified  in  Eq.  (6).  This  approach 
easily  extends  to  redundant  systems  where  J  may  not  be  square. 

To  understand  how  this  technique  handles  the  singularities 
of  the  Jacobian  J2  for  the  rotational  simulation  problem,  the 
objective  in  this  case  is  to  minimize  L  no>v  being  specified  by; 


L  =  IQ-y2«l’«,  +  X-|@l’^;  (12) 

where  X  >  0  is  selected  to  minimize  joint  velocities  if  the 
Jacobian  becomes  singular.  One  could  rewrite  Eq.  (12)  to  ana¬ 
lyze  each  term  separately.  Denote  a  residuals  term  r  1  and  an 
energy  term  r2  with  weighting  matrices  Rl  and  R2  as  follows: 

(rl)-=in-/2el>.. /?1>0  (13) 

(/■2)‘  =  X‘iei'jK-, /?2>0  (14) 

The  Damped  Least  Squares  analogy  related  to  the  above 
problem,  involves  the  selection  of  0  such  that  the  following 
relationship  is  approximately  correct: 

e=.iy2r‘n  (15) 

Thus  the  fidelity  of  the  motion  field  (denoted  by 
IQ  -JlOt'm  )  is  compromised  by  trading  off  a  cost  associated 
with  a  weighting  in  By  minimizing  Z.  of  Eq.  (12)  for 

the  rotational  motion  study,  this  yields  the  augmented  equa¬ 
tion: 

col[/2,  X/]«  =  col(Q,01  (16) 

which  requires  solving: 

{J2^J2  +  \^D8  =  nTQ  (17) 

or 

e  =  (72^/1  (18) 

as  the  inverse  kinematics  solution  for  the  rotational  simulation 
problem  which  is  similar  in  form  to  the  linear  acceleration 
case. 

This  procedure  will  now  be  worked  with  an  example  from 
a  motion  simulation  problem  on  a  centrifuge  to  illustrate  the 
applicability  of  this  approach. 

Application  of  the  Approach  to  a  Motion  Simulation 
Problem 

The  concept  of  treating  a  very  large  motion  simulator  (cen- 
trigue)  as  a  “large  robot”  leads  to  some  interesting  new  per¬ 
spectives  on  how  to  dynamically  control  such  a  device.  In 
Repperger  (1990)  it  was  shown  that  work  volume  (or  perform¬ 
ance  volume)  plots  of  such  a  device  can  be  obtained  in  linear 
acceleration  space,  rotational  velocity,  and  rotational  accel¬ 
eration  space.  This  can  be  related  to  an  acceleration  radius 
concept  (Graettinger  and  Krogh,  198R>  in  which  a  spherical 
volume  figure  is  obtained  which  represents  the  (three  dimen¬ 
sional)  linear  accelerations  attainable  at  the  end  effector.  It  is 
noted  that  this  type  of  volume  is  restricted  to  a  spherical  figure; 
the  performance  volume  in  (R'^  ,jT>erger,  1990),  however,  is  not 
restricted  to  any  particular  geometric  shape.  Since  performance 
volumes  of  this  type  can  te  determined,  the  next  step  in  un¬ 
derstanding  belter  ways  to  control  this  simulator,  if  a  motion 
trajectory  is  specified  apriori,  involves  the  solution  of  the  in¬ 
verse  kinematics  problem.  To  accomplish  this  task,  the  deter¬ 
mination  of  the  joint  motions  is  required  to  emulate  a  specified 
path  trajectors-  ot  some  endpoint  within  the  centrifuge.  The 
last  step  for  riie  control  engineer  to  regulate  the  motion  of  this 
simulator  would  be  to  solve  the  inverse  dynamics  problem. 
This  in'  oives  the  calculation  of  the  required  torques  (or  com¬ 
mand  inputs)  at  the  actuators  to  produce  the  complex  accel¬ 
eration  field  desired  at  a  point  of  interest  within  the  centrifuge. 
In  this  study  the  required  torques  arc  bounded  by  the  variable 
C  in  Eq.  (7)  (and  thus  the  inverse  dynamics  problem  is  solved 
as  a  consequence  of  this  procedure)  using  the  Damped  Least 
Squares  procedure.  This  method  implicitly  utilizes  the  con¬ 
straint  relationship  in  Eq.  (7).  First  the  forward  kinematics 
problem  must  be  formulated  to  determine  the  Jacobian  J2  and 
to  observe  if  the  motion  simulation  problem  is  well  posed. 

Problem  1  -  The  Forward  Kinematics  Problem 

Figure  I  illustrates  a  diagram  describing  the  DES  (Dynamic 
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Fi9. 1  The  OES  (Dynamic  Environment  Simulator) 
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Fig.  2  The  body  axis  coordinate  system  used  lor  human  studies 

Environment  Simulator).  This  device  is  a  19  foot  radius  cen¬ 
trifuge  with  the  principal  axis  of  rotation  (^i)  as  shown  in  the 
figure  with  an  angular  velocity  In  general  the  arm  rotates 
clockwise  when  viewed  from  the  top  as  indicated  in  the  figure. 
At  the  end  of  the  arm  is  the  fork  which  rotates,  using  its  own 
control  system,  about  the  pitch  axis  (62).  The  cab  also  rotates 
independently  about  the  same  center  of  rotation  as  the  fork 
but  in  the  roll  axis  as  illustrated  in  Fig.  1 .  The  human  subject 
is  positioned  in  the  cab  as  indicated. 

Within  the  cab,  at  the  seat,  is  illustrated  the  body  centered 
coordinate  system  used  in  acceleration  research  for  human 
subjects  who  ride  centrifuges  (Fig.  2).  This  includes  a  yaw 
rotation  of  the  seat  about  the  body  centered  { -  z)  axis  (coin¬ 
cident  with  the  vertical  axis  of  the  cab),  and  a  seat  back  angle 
6,.  The  objective  is  to  calculate  the  accelerometer  readings  at 
point  Pj  in  Fig.  3  induced  by  the  simulator  motion  only.  In¬ 
itially  the  force  of  gravity  is  not  considered,  however,  when 
the  forward  kinematic  chain  is  determined,  the  directional 
cosines  of  the  transformation  matrix  will  indicate  the  contri¬ 
bution  of  the  gravity  vector  to  the  readings  of  the  accelero¬ 
meters  positioned  on  the  seat  at  point  P5. 

In  Fig.  3  at  the  point  Pj  of  the  body  centered  coordinate 
system  located  on  the  seat  displayed  in  Fig.  2,  it  is  first  desired 
to  calculate  the  linear  acceleration  experienced  with  respect  to 
all  the  variables  W|,  «],  LI,  L2,  $2,  $2,  ©2,  9j,  Bj,  h.  Lx, 
Ly,  Lz,  Lt,  B„  and  It  is  noted  that  the  variables  LI,  L2, 
(defined  in  Fig.  !),  Lx,  Ly,  Lz,  Lt,  B„  and  B,  (defined  in  Fig. 
3)  are  fixed  quantities  which  never  change  once  the  simulator 
motion  begins.  Table  1  lists  the  actual  values  of  these  fixed 
variables  as  they  apply  to  the  motion  simulator  in  Fig.  I  with 
the  maximum  and  minimum  values  of  the  rotational  variables 
listed. 

To  develop  the  forward  kinematic  relationships  due  only  to 
simulator  motion,  it  is  required  to  use  6  coordinate  systems. 


Fig.  3(W 

Fig.  3  The  coordinate  systems  used  lor  the  Denavit-Hartenberg  lor. 
mutation  ol  the  lorward  kinematics  equations 


Let  the  point  Pg  be  at  the  center  of  the  coordinate  system  (A’o, 
To,  Zo)  and  attached  to  the  stationary  frame  of  the  centrifuge 
(the  rigid  body  called  the  fixed  base,  (cf.  Fig.  (3)),  the  point 
Pi  is  attached  to  the  arm  which  is  also  illustrated  in  Fig.  3  as 
being  the  origin  of  the  coordinate  system  (X,,  Ki,  Zi)  (the 
second  rigid  body).  The  third  coordinate  system  (A'j,  Tj-  ^2) 
at  point  P2  is  attached  to  the  fork  of  the  centrifuge  at  its  center 
of  rotation.  The  fourth  coordinate  system  ((Xj,  Tj,  Zj)  at  point 
Pj)  is  attached  to  the  cab  at  the  same  center  of  rotation  of  the 
fork/cab.  The  point  P4  (with  coordinate  system  (X^,  K».  Z4)) 
is  a  translation  of  ( -  Lx,  -  Ly,  Lz)  units  to  the  bottom  the 
seat  (at  floor  level).  The  point  Pj  is  a  translation  of  Lt  units 
up  the  seat  (after  the  seat  is  rotated  -  6,  degrees  about  -  Z4') 
to  the  point  where  the  accelerometers  are  mounted.  Table  2 
illustrates  the  D -//  parameters  for  the  centrifuge  simulator. 

To  determine  the  Jacobian  Jl  (for  the  linear  stress  problem 
previously  introduced),  the  last  coordinate  system  (cf.  Fig.  3) 
of  interest  is  written  (cf.  Craig,  1986): 

coUJfl?,  Yll,  Z171=  ro”col[,i(0,  1^,  ZOl  (19) 
and  can  be  expressed: 
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Table  1 

The  fixed  quantities  of  the  centrifuge  simulator 

Fixed 

Variable 

Dimensions 

Rotational 

Variable 

Ma-x  .Min 

Dimensions 

II 

10.4  ft 

0, 

6.28  0 

radians 

LI 

19.0  ft 

9, 

5.86  -5.86 

radians/  second 

Lx 

3.3  ft 

Q, 

0.25  -0.25 

radians  ^s/s 

Lv 

0.1  ft 

82 

6.28  0 

radians 

Lz 

0.1  ft 

O2 

3,14  -3.1-' 

radians /second 

Lt 

2.3  ft 

a. 

2.5  -2.5 

radians  s/s 

6. 

1.57  radians 

9, 

6.28  0 

radians 

.53  radians 

9i 

3.14  -3.14 

radians/ second 

n. 

5.0  -5.0 

radians/s's 

Table  2 
fuge 

The  Denavit-Hartenberg  parameters  for  the  centri- 

A.ngle  From  Distance  Along  Distance  Along  From 

Link  c 

/-!  to  zi 

zi- 1  to  zi 

Xi-\  to  Xt  Xi-l  to  Xi 

about  Xi-I 

along  A'- 1 

along  Zi 

about  Zi 

i 

-Qi-\ 

=  Ai-l 

-di 

=  01 

1 

0 

0 

0 

-81 

2 

0 

L2 

0 

0 

T, 

0 

0 

0 

tc/2 

4 

-*•  ir/2 

0 

0 

0 

5 

0 

0 

0 

82 

6 

0 

0 

0 

+  t/2 

7 

—  t/2 

0 

0 

0 

8 

0 

0 

0 

83 

9 

0 

-  Lx 

*Lz 

0 

10 

0 

0 

0 

-  t/2 

11 

0 

Ly 

0 

0 

12 

+  ir/2 

0 

0 

0 

13 

0 

0 

0 

14 

-t/2 

0 

0 

0 

15 

0 

0 

0 

-8t 

16 

*T/l 

0 

0 

0 

17 

0 

0 

-Lt 

0 

(R  Px\\ 


(20) 


where  /?  is  the  3  x  3  rotation  matrix.  The  Jacobian  can  be 
expressed: 


jdPx\/dB\, 

hf 


1. 


BPyl/dBl, 


Pzl/dBl, 


dPxl/dB2, 
BPyi/dei, 
dPz  1/382, 


dPxi/ees: 

BPyi/dei, 

dPzl/de3,^ 


(21) 


The  elements  of  the  matrices  in  Eqs.  (20)  and  (21)  are  not 
written  out  in  detail  due  to  their  complexity.  The  derivation, 
however,  is  straightforward  as  described  using  Eqs.  (19)-(21). 


Accelerometer  Readings  Due  Only  to  Gravity 

Once  the  transformation  matrix  Tg’’  is  derived  from  the 
Denavit-Hartenberg  parameters  describing  the  centrifuge  sim¬ 
ulator  in  Table  2,  the  linear  accelerations  due  only  to  gravity 
can  be  determined  from  the  directional  cosines.  To  determine 
these  readings,  one  must  find  the  component  of  gravity  in  the 
final  coordinate  system  with  unit  vectors  (if,  k^).  Due  to 
the  fact  that  the  To' '  has  an  orthogonal  inverse  matrix  Tn°,  it 
can  easily  be  shown  that  the  gravity  vector  produces  an  ac¬ 
celeration  in  the  last  coordinate  frame  with  the  following  com¬ 


ponents: 

;5COmponent=  -  ri7°(3,  1)  (22) 

=  -s2  sy  ct  +  c2  c3  st  -  c2  s3  cy  ct  (23) 

y5COmponent=  -  rn®(3,  2)  (24) 

=  -s2  cy-c2  s3  sy  (25) 


^jcomponent  =  -  3) 

=  ~sy  sr  s2  +  c2  c3  ct  c2  s3  cy  st  (27) 

where  si  and  ci  are  notation  for  sine  and  cosine.  Equations 
(22)-(27)  describe  the  readings  of  the  accelerometers  due  only 
to  the  1  G  gravity  vector  at  the  point  P^  within  the  centrifuge 
which  act  on  the  pilot  as  he  sits  in  the  seat  where  the  endpoint 
was  previously  defined.  The  next  s’  p  is  to  calculate  the  linear 
and  rotational  accelerations.  These  accelerations  produce  forces 
and  torques  which  act  on  a  pilot  which  stress  him  in  a  phys¬ 
iological  sense. 

The  derivation  of  the  linear  accelerations  can  now  be  ac¬ 
complished  by  first  expressing  a  vector  from  (A'u,  Yg,  Zj) 
to  (Xf,  Ys,  Zj)  and  twice  differentiating  this  vector.  .Note  that 
if  P;  is  expressed  in  terms  of  the  unit  vectors  of  the  inertial 
frame  (ig.  jg,  kg),  then  by  differentiating  this  expression,  the 
time  derivatives  of  the  unit  vectors  (lo,  jg,  kg)  vanish.  To  de¬ 
termine  Pi,  it  is  composed  first  of  a  vector  P,  from  iX,,.  Yg, 
Zg)  to  (Xi.  Yi,  Zj),  which  can  be  written  as: 

Pj  =  i.2/,  (28) 

and  a  vector: 

P4  =  -  Lx  1}  +  Lz  k;  -  Ly  j-.  -h  P,  (29) 

which  describes  the  translation  from  the  points  Pi  to  P4.  Finaily 
a  vector  defined  by: 

Pi--^Pi  +  P„f-Pi-Lt  ki  (30) 

describes  the  last  position  change  to  the  final  coordinate  sys¬ 
tem,  With  some  work  (and  following  a  procedure  similar  to 
expressing  (15,75,  k^)  in  terms  of  the  unit  vectors  of  the  inertial 
system  {tg,  jg,  icg)),  the  following  relationships  can  be  shown 
to  be  true: 

/i  =c!  lo-sl  7o  (31) 

i)  =  (  -5ls2c3  -  els3)  to  +  (  ~  c’l'S2c3  -t-  5ls3)yo  -  ( c2c3)  kg  (32) 

yj=  (sl52s3  -clc3)  ig+  (clj2a3  -“5103)70 (  -  cZsJ)  kg  (33) 

kx=  (  -5lc2)  ig*  (  -cIc2)7o  +  (-52)  kg  (34) 

Ar5  =  (r,7°(l,3))fo+  (7-r®(2.3))7o+  (rr“(3.3))ko  (35) 

where  the  Tp®  components  are  defined  as  the  orthogonal  in¬ 
verse  of  the  matrix  in  Eq.  (20),  Thus  expressing  P5  in  terms 
of  the  inertial  frame  yields: 

P5  =  {AO{t))io+  (B0(()  )7o+  (CO(O)^o  (36) 

where  A0(t),  ff0(0,  and  C0(0  are  written  out  in  some  detail  in 
(Repperger,  1990).  This  expression  can  now  be  differentiated 
twice  with  respect  to  time  yielding: 

Ps=  (/ll(/))/o+  (Bl(0)jo+  (Cl(O)ko  (37) 

where  the  expressions  for  <4 1(0,  P1(0,  and  C1(0  are  given  (for 
completeness)  in  (Repperger,  1990).  The  next  step  in  this  pro¬ 
cedure  is  to  develop  the  rotational  relationships  and  the  Ja¬ 
cobian  J2  of  the  centrifuge  simulator. 

Derivation  of  the  Rotational  Relationships 

Determination  of  the  Jacobian  J2  for  the  Centrifuge:  To 

apply  these  techniques  for  the  supermaneuverability  flight  tra¬ 
jectory  to  a  centrifuge  simulator,  first  it  is  necessary  to  deter¬ 
mine  both  the  Jacobian  term  J2  and  define  the  rotational 
trajectory  of  interest.  At  the  point  Pj  at  the  endpoint  of  the 
end  effector  (the  point  on  the  seat  within  the  centrifuge),  the 
total  angular  velocity  of  the  third  fixed  body  (the  cab)  can  be 
written  in  vector  form  as; 

-  81  k0+B2  k2+B3  k3  (38) 
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Fig.  4  The  “HerbsCs”  maneuver— e  supermaneuver 


where  i,  j,  k  represent  unit  vectors  in  the  coordinate  systems 
defined  in  Fig.  3.  It  is  noted  that  the  angular  velocity  of  the 
human  subject  within  the  cab  is  identical  to  the  angular  velocity 
of  the  cab  (because  the  seat  is  fixed  throughout  the  run).  It  is 
necessary  to  express  kl  and  kl  in  terms  of  the  inertial  coor¬ 
dinates  (iO,jQ,  kO).  From  the  relationships  derived  in  Table  2, 
it  is  seen  that: 

^2  =  (cl)/0-(sl)/0  (39) 

and  kl  can  be  described  via; 

/!r3  =  (-j/c2)/0-t-(-cl  c2)^  +  (-s2)«)  (40) 

Now  substituting  (39)  and  (40)  into  (38)  yields  for  the  fi  of  the 
third  rigid  body  to  be  expressed  in  terms  of  the  inertial  frame 
as; 

fl  =  /12(t)i0  +  B2(t)j0  +  C2(0«>  ('tl) 

where  for  completeness,  it  is  seen: 

/12(f)  =  91c\  -  93  si  c2  (42) 

B2{t)  =  -  92sl  -  h  cl  c2  (43) 

C2(/)= -ei-93  52  (44) 

To  determine  the  Q  in  the  last  coordinate  system  at  the  subject, 
Eq.  (4!)  must  be  multiplied  by  7g'’  from  (Repperger,  1990, 
1991)  to  yield; 

n  =  /I3(/)/0  + 53(f)/) +  C3(f)«)  (45) 

where  it  can  be  shown  that; 

/43(f)  =  hcl  -  92^51  -  '9yslc2  -  AjS iclc2  +  « (46) 


Fig.  5  An  approximation  for  the  roll,  pitch,  and  yaw  lime  histories  with 
respect  to  the  aircraft’s  axes  when  emulating  the  Herbst's  maneuver 

Application  to  the  First  Supernianeuver — The  Herbst’s 
Maneuver 

This  is  the  first  classical  supermaneuver  that  was  designed 
from  the  original  concept  of  Herbst  (Herbst,  1972).  However 
this  method,  it  turns  out,  will  be  applied  to  a  trajectory  which 
produces  a  singular  Jacobian  matrix.  To  specify  0(f)  for  the 
supermaneuverable  trajectory  called  the  “Herbst’s  Maneuver” 
(Herbst,  1972,  1980,  1983)  as  illustrated  in  Fig.  4,  and  in  terms 
of  body  axis  rotations  of  the  airplane  in  Fig.  (5),  the  following 
approximations  are  made  for  the  roll,  pitch,  and  yaw  rotations 
(f  is  in  units  of  seconds  and  the  angles  are  in  radians); 

Roll  rotation  =  @x(f)  =  ,16  x-  sin  l(x/5)f]  (53) 

Yaw  rotation  =  9z(f)  =  (T/180)f^  for  0<f<5  (54) 

=  .436.' for  5<f<7  (55) 

=  T(l-e-''-‘  '’)  for  7  <  f<  10  (56) 

Pitch  rotation  =  5y(f)  =  x(. 22 -.19  cos  [(T/5)fJ)  (57) 

This  yields  for  rotational  velocity  space: 

Qx(t)  =  .03  cosI(t/5)  f)  (58) 

Oj'(f)  =  .38  sin(x/5)  f  (59) 

Qz(f)  =  T/60  f^  forO  <  f<5  (60) 

=  .436  for  5  <  f  <  7  (61) 

=  Tt-''-‘'‘’for7  <  f  <  10  (62) 

The  matrix  J2  can  easily  be  obtained  for  the  centrifuge  from 
the  relationshi.p  (49)  by  writing  Q  in  terms  of 
61,  62,  and  63.  The  terms  of  J2  are  listed  here  for  complete¬ 
ness: 


B3(0=  -  62SI  -  626,c1  -e3Clc2-)-  0je,slc2+  636JCI52  (47)  fl.  The  Jacobian  of  Ibe  Centrifnie 


C3(f)= -6,-63^2-6362  c2  (48) 

But  recall  from  the  definition  of  the  J2  Jacobian  term,  the 
relationship: 

Q  =  J29  (49) 

translates  into: 

6=72-'n  (50) 

when  J2  is  invertible.  For  simplicity  to  integrate  Eq.  (50)  for¬ 
ward  in  time,  a  first-order  Euler  method  will  be  used  of  the 
form: 

(51) 

or 

6**,  =  6*  +  AfI/2-']a(f)  (52) 

where  is  specified  apriori. 


/  ' 
'  sy  ct,  -St  cl  +  cy  ct  si ,  si  c2  st  +  cl  c2  cy  cf +s2  sy  ct 

-cy,  si  sy,  clc2sy-cys2 

^^-sy  St,  -cl  cf-sl  cy  st,  si  c2  cf-cl  c2  cy  sf-s2  sy  st 


The  objective  is  to  determine  an  inverse  kinematics  solution 
similar  to  Eq.  (6).  A  specific  seat  yaw  rotation  and  seat  back 
angle  have  been  assumed  (9y  =  0  degrees,  6f  =  90  degrees). 
The  Jacobian  in  this  case  represents  parameters  of  the  centri¬ 
fuge  simulator  and  is  specified  via: 


0,  -cl,  si  c2 
/2=  -1,  0,  -s2 

0,  -si,  -cl  c2 
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Fig.  6  Hsrbsl’s  maneuver  simulation 


-•-Om*ga’  Exoc^  -A- App'ox.  ; 

Fig.  7  Herbst’s  maneuver  simulation  using  damped  least  squares 

The  X(t)  vector  in  Eq.  (1)  corresponding  to  the  “Herbst’s 
Manuever”  is  described  in  Eqs.  (i8)-(62)  for  t « [0,  10].  It  can 
be  shown  (Fig.  6)  that  the  determinant  of  J2  vanishes  at  f  =  5.5. 
Applying  a  Frobenious  norm  at  r=  5.5  when  J2  becomes  rank 
deficient,  it  is  determined  that  II.YI1^  =  14.2.  The  limits  on 
joint  velocities  and  torques  in  their  appropriate  dimensions 
requires  C  <  1025.,  thus  resulting  in  (at  t  =  5.5),  ll/2iVr  = 
3.01,  even  though  J2  is  rank  deficient.  Thus  C  in  Eq.  (7) 
constrains  d  such  that: 

8  1025.0  (65) 

in  the  appropriate  joint  velocity /torque  dimensional  units  for 
this  problem.  It  can  easily  be  shown  that  if  the  selection  is 
made; 

X  =  5  (66) 

is  an  appropriate  choice  of  this  variable.  In  a  least  squares 
sense  this  will  provide: 

I  1025.0  (67) 

Figure  7  illustrates  this  implementation  for  \  =  5  with  the 
resulting  decrease  in  the  joint  velocities  constrained  via  8  6  9^ 
<  1025.0.  Notice  also  the  compromise  in  the  fidelity  of  the 
motion  field  simulation  which  has  to  be  traded  off  in  order 
to  limit  the  joint  velocities  and  torques.  Figure  6  shows  an 
attempted  solution  without  using  damped  least  squares  and 
the  apparent  divergence  of  the  computed  inverse  kinematics 
algorithm  at  f  =  5.5. 

Application  to  the  Second  Supermaneuver— The  “Cobra 
Maneuver” 

The  second  trajectory  selected  is  termed  the  “Cobra  Ma¬ 
neuver.”  This  has  been  demonstrated  by  actual  aircraft  in  flight 
scenarios  {Aviation  Week  and  Space  Technology,  1990).  For 
this  particular  supermaneuver,  the  assumption  is  made  for  the 


40 - — - — - — - — _9 


-120 - - - - - • - - - 4 

0125456789  10 

Tima  lo  Socoods 


FoFK  -A-  Arm  CoO 


Fig.  8  Cobra  simulation  on  centrifuge 


centrifuge  simulator  that  dy  =  t/2,  dt  =  0,  and  t  e  (0,  10] 
seconds.  The  angles  are  in  radians.  The  Jacobian  now  becomes 
(for  this  seat  configuration); 

1,  0,  s2  \ 

0,  jl,  clc2  (68; 

0,  -cl,  slc2i 

The  Cobra  Maneuver  was  simulated  using  the  following  an¬ 


gular  rotations: 

Roll  rotation  =  e3(r)  =  (t/180)(30-  2.50  (69) 

Yaw  rotation  =  0.0  (70) 

Pitch  rotation=  -02(/)=  (ir/3)[l  -cos(Tf/5))  (71) 

011(0  =  8. 8 -.66  t  in  RPM  for  t  e  [0,  5]  (72) 

0)1(0  =  5.5 +  .66  [r-5J  for  r  e  [5,  10]  (73) 


It  is  noted  that  the  determinant  of  J2  for  the  Cobra  maneuver 
remained  close  to  the  value  of  1.0  during  this  simulation  and 
no  possibility  of  a  singular  Jacobian  occurred.  The  expression 
(72)  for  0)1(0  is  used  to  demonstrate  the  slowing  down  of  the 
aircraft  as  it  assumes  its  unusual  attitude.  In  the  centrifuge, 
the  cab  has  to  vector  out  to  compensate  for  gravity  which  is 
reflected  by  the  &3{f)  term  to  compensate  for  the  arm  slowing 
down.  Figure  8  illustrates  this  numerical  simulation. 

Application  to  the  Third  Supermaneuver — The  High  Yo- 
Yo  Turn  For  BFM 

Of  the  43  supermaneuvers  and  38  variations  in  (Foltyne  et 
al.,  1987),  the  third  maneuver  selected  was  the  High  Yo-Yo 
Turn  for  Basic  Fighter  Maneuver  (BFM).  This  was  chosen 
because  this  supermaneuver  precludes  overshoot  and  develops 
control  closure  (maneuver  16  of  the  above  reference).  When 
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Fig.  10  Centriluge  responses  to  the  joint  commands 

an  aircraft  flies  this  trajectory,  it  maintains  controllability  dur¬ 
ing  the  combat  engagement.  Figure  9  illustrates  this  maneuver 
as  the  agile  aircraft  flies  this  trajectory  versus  his  opponent. 
For  the  centrifuge  simulation,  this  is  strictly  a  yaw  and  pitch 
motion  (assuming  the  same  seat  configuration  as  with  the 
Herbst’s  maneuver). 

This  motion  field  simulation  turns  out  to  be  easier  to  im¬ 
plement  than  the  Herbst’s  maneuver  due  to  the  longer  time 
required  to  complete  the  yaw  motion.  Mathematically  this 
maneuver  can  be  characterized  by  the  following  yaw  motion: 

(?:(/)  =  l9-(/)  =- (ir/3)(//10)  (74) 

and  the  following  pitch  motion; 

-@3(f)  =  9,(r)  =  (5r./180)[20-(-50  sin  [2rr //IS]  (75) 

The  determinant  of  the  Jacobian  is  better  behaved  as  compared 
to  the  Herbst’s  maneuver  since  it  remains  greater  than  zero 
throughout  the  motion  trajectory.  Figure  10  illustrates  the 
commanded  and  actual  responses  of  joints  2  and  3  for  this 
supermaneuver.  Finally,  some  comments  need  to  be  made  on 
the  inverse  dynamics  problems  associated  with  these  simula¬ 
tions. 

The  Issues  Related  to  Inverse  Dynamics 

From  observing  Fig.  10,  more  discussion  is  required  to  de¬ 
scribe  the  inverse  dynamics  problem  in  the  simulation  of  these 
supermaneuverable  trajectories.  Space  precludes  an  extensive 
discussion  in  this  area;  however,  it  is  important  to  state  the 
following  facts: 

(1)  The  arm  (wl(0)  is  a  velocity  control  system. 

(2)  The  fork  and  the  cab  and  93(f))  are  both  position 
control  systems. 

(3)  Each  control  system  was  modeled  somewhat  independ¬ 
ently  with  a  disturbance  input  to  characterize  the  off  axes 
Coriolis,  centrifugal,  and  unmodeied  gravity  terms. 

(4)  Step  responses  have  been  recorded  in  each  axes  at  various 
arm  speeds  to  obtain  linearized  models  about  specific  operating 
points. 

The  implementation  of  the  solution  to  the  inverse  dynamics 
problem  can  be  illustrated  in  maneuver  number  3  for  the  92(f) 
response  (cf.  Fig.  10).  Since  this  is  a  position  control  system. 


the  62  input  commands  arise  from  the  inverse  kinematic?  so- . 
luiions.  The  actual  92  response  of  the  simulator,  however, 
differs  from  the  92  input  commands.  A  closed-loop  position 
control  system  is  used  to  drive  the  fork  with  a  forward  loop 
gain  K.  This  closed-loop  system  includes  the  dynamics  of  the 
fork  and  a  disturbance  input  to  account  for  unmodeled  effects 
and  the  undesirable  Coriolis  terms.  The  gain  K  is  increased  to 
reduce  this  position  error  but  still  must  be  kept  small  enough 
to  maintain  closed-loop  system  stability.  There  is  an  obvious 
tradeoff  between  stability  and  performance  in  this  case  and 
the  issues  related  to  the  tradeoff  between  the  stability  of  this 
inner  loop  to  the  resulting  compromise  of  the  fidelity  of  the 
motion  field  will  be  discussed  in  future  papers  on  this  topic. 
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